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Abstract 

We discuss recent work with E.Floratos (JHEP 1004:036,2010) on Nambu Dy- 
namics of Intersecting Surfaces underlying Dissipative Chaos in R 3 . We present our 
argument for the well studied Lorenz and Rossler strange attractors. We implement 
a flow decomposition to their equations of motion. Their volume preserving part 
preserves in time a family of two intersecting surfaces, the so called Nambu Hamil- 
tonians. For dynamical systems with linear dissipative sector such as the Lorenz 
system, they are specified in terms of Intersecting Quadratic Surfaces. For the case 
of the Rossler system, with nonlinear dissipative part, they are given in terms of a 
Helicoid intersected by a Cylinder. In each case they foliate the entire phase space 
and get deformed by Dissipation , the irrotational component to their flow. It is 
given by the gradient of a surface in R 3 specified in terms of a scalar function. 
All three intersecting surfaces reproduce completely the dynamics of each strange 
attractor. 

1 Introduction 

Dissipative dynamical systems, with a low dimensional phase space, present an important 
class of simple non-linear physical systems with intrinsic complex behavior (homoclinic 
bifurcations, period doubling, onset of chaos, turbulence), which generated intense exper- 
imental, theoretical and numerical work in the last few decades [TJ [2] . 

Recently [3] we have reexamined dissipative dynamical systems with a 3-dimensional 
phase space from the perspective of Nambu-Hamiltonian Mechanics(NHM) [U [5]. The 
latter represents a generalization of Classical Hamiltonian Mechanics, mostly appropriate 
for the study of odd- dimensional phase-space volume preserving flows ( Liouville's theo- 
rem). As such and in order to make it directly applicable to the dynamics of dissipative 
systems in R 3 we must associate it to a volume preserving dynamics sector. We have done 
so by introducing a flow decomposition to their equations of motion and therefore isolate 
in their flow vector field its rotational(solenoidal) part. It is manifestly volume preserv- 
ing, non-dissipative and hence directly describable in terms of the intersecting surfaces 
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of NHM. Along with the remaing irrotational component to the flow the decomposition 
fully recovers the dissipative dynamics of the system in question, presenting itself as an 
equivalent formulation of the odd- dimensional dynamical system. 

We have applied the above recipe to the famous examples of Lorenz [6] and Rossler [7] 
chaotic attractors which represent the prototype models for the onset of turbulence. [2]. 

In sect. 2 we start off with a discussion of flow decompostion for the most elementary 
and familiar of all dissipative systems in R 2 without chaotic behavior : the dissipative 
harmonic oscillator (DHO). In this case the existence of a Hamiltonian formalism identifies 
for its "closed" nondissipative sector the harmonic oscilator (HO) with a well defined 
integrable classical and quantum evolution. 

The presence of chaotic flows for dissipative systems in R 3 is argued, by analogy, 
to nescessitate intersecting Nambu Surfaces which define rerspectively "closed" physical 
systems with integrable R 3 periodic orbits and simultaneously a well defined classical and 
quantum behaviour. 

In sect. 3 we apply the framework for the cases of the Lorenz and Rossler strange 
attractors. We isolate their non-dissipative sector parametrized by two intersecting sur- 
faces: a cylinder and a paraboid for the Lorenz attractor as well as a helicoid with a 
cylinder for the Rossler system. They account, amazingly, for the double scroll topology 
of the full " butterfly" Lorenz attractor. and the single scroll topology for the Rossler 
case. We end our presentation with Conclusions and open problems. 

2 Flow Decomposition in Dissipative Systems: R 2 
versus R 3 

Dissipation is a necessary condition for dynamical systems to exhibit chaos. Yet it is not 
suficient. This is enunciated through a powerful No-go theorem . Indeed the Poincare- 
Bendixon alows for only fixed points and limits cycles in two dimensions. Chaotic Flows 
need space to emerge. At the minimal level they emerge with dissipative systems in 
R 3 . They are typically associated with Strange Attractors such as the famous ones of 
Lorenz and Rossler. They belong to a large class of dynamical systems whose dynamics 
is governed by continuous set of 1st order ordinary differental flow equations 

x = v(x(t), t, A) (1) 

where v is a velocity field flow with A some control external parameter. Their phase space 
dynamics, depending on whether they exchange energy with their environment or not, 
can be either open- dissipative or closed (conservative- Hamiltonian) . This is reflected on 
their velocity flow field being divergenceless or not. Let us see all these issues for the 
simple case of the Dissipative Harmonic Oscillator (DHO) whose phase space equations 
of motion are well known to be for X\ — q and X2 = p: 



q = p + aq 
p = —q + /3p 



(2) 



The velocity flow is given by Vv tot = (a + (3) ^ for appropriate values a ^ —j3 ^ with 
a net outflow or inflow of energy depending on the sign of Vvtot (inflow < or outflow 
> 0). Nevertheless we split it into a nondissipative component v^d = (p, —q) with zero 
flow Vvnd = 0. Its dissipative part is given by = otq + (3p with its flow being simply 
the total flow of the HO Vvd = a + /3. The dynamics of the HO in one dimension can 
take the form of a 2nd order differential equation 

q+(a + + (1 + = (3) 

The damping factor is 7 = a + (3 and the effective natural frequancy which is dissipation 
induced is given by u e f f = 1 + a (3 . 

The existence of the Harmonic Oscillator Hamiltonian H — \{p 2 + q 2 ) along with a 
Dissipation function D = \{aq 2 + f3p 2 ) reproduce the eqs. of motion of (2.2) in a more 
compact form 

•J-' = c'HOjH-0,1)) (4) 

where i,j = 1,2 and e 12 = — e 21 = 1 as usual. The non-dissipative sector of the HO is 
given by 



dH 

4 = -TT~ = V 

op 

V = = -q (5) 

Firstly it identifies a circular periodic orbit, the famous HO as the integrable progenitor 
of the transient DHO . This is a well defined closed physical system with a most familiar 
quantum behavior, the Quantum Harmonic Oscillator. 

Secondly the linear dissipation of the DHO a transient spirals it inwards to a fixed point 
of zero energy. The harmonic oscillator localizes in effect such a dissipative evolution . 

Lastly the rate at which the energy of the harmonic oscillator loses its energy as 
it spirals in depends on its damping strength. It is the Quality factor which is easily 
computed to be 

7 oc + p 

All in all Flow decomposition for the transient DHO has three immediate implications for 
the DHO: 1. Existence of Integrable Harmonic oscillator associated with a periodic orbit 



2. Localization of the DHO evolution and asymptotic fixed point state 3. Quality factor 
for DHO as a measure of its energy loss rate and in effect of the damping strenth. 

We will proceed now to examine whether it is feasible to implement this methodology 
for the case of Strange attractors which possess chaotic flows in R 3 . 

Nambu-Hamiltonian mechanics is a specific generalization of classical Hamiltonian 
mechanics, where the invariance group of canonical symplectic transformations of the 
Hamiltonian evolution equations in 2n dimensional phase space is extended to the more 
general volume preserving transformation group SDif f(M) with an arbitrary phase space 
manifold M of any dimension d = dim(.M) . 

In [3] we work with the case of a three dimensional flat phase space manifold. Never- 
theless our results could be generalized to curved manifolds of any dimension [S]. 

Nambu-Hamiltonian mechanics of a particular dynamical system in R 3 is defined once 
two scalar functions Hi £ C°°(R 3 ),i = 1,2, the generalized Hamiltonians [H [5] are pro- 
vided. The evolution equations are: 

x i = {x\H u H 2 } z = 1,2,3 (7) 

where the Nambu 3-bracket, a generalization of Poisson bracket in Hamiltonian mechanics, 
is defined as 

{f,g,h} = & k &f&gdFh i,j,k = 1,2,3 V/, g, h £ C°°(R 3 ) (8) 
Any local coordinate transformation 

x* -> y l = y\x) z = l,2,3 (9) 

which preserves the volume of phase space 

d6t (Iz^) = 1 V x = (x l ,x 2 ,x 3 ) e R 3 (10) 

leaves invariant the 3-bracket and therefor it is a symmetry of Nambu-Mechanics. Except 
for the linearity and antisymmetry of the bracket with respect to all of its arguments it 
also satisfies an important identity, the so called "Fundamental identity" [FI] [21 E]. 
The evolution eq.(2.7) has a flow vector field: 

v \x) = e ijk d j H 1 d k H 2 i,j,k = 1,2,3 (11) 

which is volume preserving 

<9¥ = (12) 

The reverse is also true. We will name the phase-space volume preserving flows "Non- 
dissipative" while the non-conserving ones " Dissipative" (d l v l (x) > 0(< 0)). 

The flow equations of a general dissipative system may take the general vector form 



'x = VH X x VH 2 + VD (13) 
We notice that given a pair of functions if i , H 2 such that 

v ND = V#i x VH 2 (14) 

any transformation of Hi, H 2 

Hi^H[{H x ,H 2 ) 2 = 1,2 (15) 

with unit Jacobian 
gives also 

v ND = VH'i x VH' 2 (17) 

In ref.jB] we reduced the evolution equation of the form (2.32) Nambu Mechanics in 
Hamiltonian-Poisson form as follows: 

x* = {x\H u } H2 (18) 

where the induced Poisson bracket 

{f,g} H2 = <■ >>-<)■ .(<>><,<)'■■ n 2 (19) 

satisfies all the required properties like linearity, antisymmetry and the Jacobi identity. 

In the following sections we are going to present a detailed investigation of the Lorenz 
and Rossler attractors from the point of view of Dissipative Nambu-Hamiltonian Dynam- 
ics. 



3 The Lorenz Attractor from Dissipative Dynamics 
of Intersecting Quadratic Surfaces . 

The Lorenz model was invented as a three Fourier mode truncation of the basic eqs. for 
heat convection in fluids in Reyleigh-Benard type of experiments [12] The time evolution 
eqns. in the space of three Fourier modes x, y, z which we identify as phase-space are : 

x = a(y — x) 
y = x{r — z) — y 

z = xy — bz (20) 



where a is the Prandtl number, r is the relative Reynolds number and b the geometric 
aspect ratio . 



The standard values for a, b are a = 10, b = | with r taking values in 1 < r < oo. 
There are dramatic changes of the system as r passes through various critical values 
which follow the change of stability character of the three critical points of the system Pi 
: x = y = z = 0, P± : x = y = ±^Jb(r — 1), z = r — 1 

Lorenz discovered the non-periodic deterministic chaotic orbit for the value r = 28, 
which is today identified as a Strange Attractor with a Hausdorff dimension of (d = 
2.06) [13J. Standard reference for an exhaustive numerical investigation of the Lorenz 
system is the book by Sparrow [H] . 

There have been various attempts made to localize the Lorenz attractor, by convex 
surfaces, in order to get information about Hausdorff dimensions [13] and other charac- 
teristics [T5] . 

We will proceed to exhibit Localization of the full Lorenz attractor from the Nambu 
surfaces which we will determine shortly. At the Quantum level the existence of an 
attracting ellipsoid in a matrix model formulation of the Lorenz system is manifest [3] as 
one gets attracting ellipsoids in higher dimensional phase spaces. 

We now proceed to describe the Lorenz system in the framework of section 2. The 
flow vector field v is analyzed into its dissipative and non-dissipative parts as follows: 

v D = {-ax,-y,-bz) = VD (21) 

with the "Dissipation" function 

D — — - (ax 2 + y 2 + bz 2 ) (22) 

and 

v N d = {(7y,x(r - z),xy) = (0, y, z - r) X (-x, 0, a) (23) 
The two Hamiltonians or Clebsch-Monge potentials if i, if 2 are determined by 

V#i x Vi? 2 = vnd (24) 

or equivalent ly 



Hi = \[y 2 + (z - r) 2 } (25) 



and 



x 2 



H 2 = az - — (26) 
The Lorenz system (20) can thus be written in the equivalent form in terms of r — (x, y, z): 

f = VFi x Vii 2 + VD (27) 
In the Non-Dissipative part(ND) of the dynamical system 

'r = VH X x VH 2 (28) 



the Hamiltonians Hi, H 2 are conserved and their intersection defines the ND orbit. More- 
over if we get the reduced Poisson structure (sect. 2) from H 2 we obtain the 2-dim phase 
space £ 2 to be the family of parabolic cylinders with symmetry axis the y-axis: 

H 2 = constant = H 2 {f Q ) (29) 

X 2 is thus given by 

z = z + ^-^ (30) 

2,(7 

with x ,z the initial condition for x,z. The induced Poisson algebra (rel.19) is given by 

{x, y} H2 = d z H 2 = a 
{y, z}h 2 = d x H 2 = -x 
{z,x} H2 = 

The dynamics on the 2d-phase space £ 2 is given by Hi 

x = {x,Hi} H2 

y = {y, Hi}h 2 

z = {z,Hi} H2 (32) 

and Hi is an anharmonic oscillator Hamiltonian with (x/a,y) conjugate canonical vari- 
ables. Using rel. (25-26) we get on £ 2 : 

Hi = ± [V 2 + \ (x 2 - « 2 ) 2 ] (33) 

with 

a 2 = x 2 - 2a(z Q - r) = -2H 2 + 2ar (34) 

where \ plays the role of the mass. Depending on the initial conditions we may have a 
sing le well (a 2 < 0, H 2 > ar) or a double well potential ( a 2 > 0, H 2 < or ) respectively. 
The trajectories , the intersections of the two cylinders , Hi and H 2 with orthogonal 
symmetry axes (x,y) may either have one lobe left/right or may be running from the 
right to the left lobe. This is reminiscent of the topology structure of the orbits of the 
Lorenz chaotic attractor. 

The full Lorenz system does not conserve Hi,H 2 and there is a random motion of 
the two surfaces against each other. Their intersection is time varying. In effect at every 
moment the system jumps from periodic to periodic orbit of the non-dissipative sector. 
Moreover the motion of the non-dissipative system around the two lobes, either left or 
right, can now jump from time to time from one lobe to the other. 



(31) 



4 The Rossler Attractor from Dissipative Dynamics 
of a Cylinder Intersecting with an Helicoid 



Rossler introduced a simpler than Lorenz's nonlinear ODE system with a 3d- phase space, 
in order to study in more detail the characteristics of chaos, which is motivated by simple 
chemical reactions [TJ. 

The Rossler system is given by the evolution eqns: 



with parameters a, 6, c usually taking standard values a = b = 0.2, c = 5.1 or a = b = 
0.1, c = 14 for the appearance of the chaotic attractor. 

We turn now to the study of the Rossler system as a Dissipative Nambu-Hamiltonian 
dynamical system. 

The key difference with the Lorenz attractor is that the dynamics of the system is 
simpler. Chaos appears as random jumps outwards and inwards the single lob attractor. 

In order to get the three scalars, the two generalized Hamiltonians if i, if 2 which are 
conserved and characterize the non-dissipative part and D the dissipation term : 



we checked after some guess work, that we must subtract and add a new term in the first 
equation. Indeed we find for the two parts 



x = —y — z 
y = x + ay 
z = b + z(x — c) 



(35) 



r = 



V#i x Vif 2 + V£> 



(36) 



VND 



-y - z - tt 




(37) 



Vd = (—,ay,z(x-c)) 

satisfying accordingly V • v^d =0 and V x v D = 
We must determine Hi, H 2 and D such that 



(38) 



V.ffi x VH 2 = vnd 



(39) 



and 



VD = v D 



(40) 



For D we find easily 

D = X - [ay 2 + {x- c)z 2 } 



(41) 



To get H 1 ,H 2 we must integrate first the Non-dissipative system: 

z 2 

x = —y — z 

y 2 

y = x 

z = b (42) 

The general solution is : 

z 2 

x(t) = -6(1 + z(t)) + (x + b(l + z )) cos(t) - ( Vo + z + -± - b 2 ) sin (t) (43) 
z 2 (t) z 2 

y (t) = b 2 - z(t) - + ( Vo + z + ^ - b 2 ) cos(t) + (x + b + bz ) sin (t) (44) 

z (t) = bt + z (45) 
To uncover Hi, H 2 we introduce the complex variable 

w(t) = wi(t) + iw 2 (t) (46) 

with 



We obtain 
with 



wi(t) = x(t) + b(l + z(t)) 

w 2 (t) = y(t) + z(t) + Z ^-b 2 (47) 



w (t) = Wo-e 11 (48) 



w = w(t = 0) (49) 
We see that there are two constants of motion, the first one being : 

I w 0) I = I w o I (50) 

and we define correspondingly, 

Hi = \ \w{t)\ 2 = l -{x + b{l + z)) 2 + \{y + z + ^-b 2 ) 2 (51) 

The second integral of motion is obtained through the phase 

w(t) =\w Q \ -e llp ° -e lt =\w \ e Mt) (52) 

or from (4.14) 

¥(t) ~ I = J ( 53 ) 



and we define appropriately the second constant surface H 2 



H 2 



b arctg 



y + z + ^ - b 2 
1 + 6(1 + 2) 



— z 



(54) 



We easily check that rel.(41) is satisfied: 



V#i x WH 2 



(-y ~ z 




(55) 



The family of surfaces Hi and H 2 are a quadratic deformation of a cylinder and 
respectively a quadratic deformation of a right helicoid. Their intersection is the trajectory 
(43-44). 

5 Conclusions-Open Problems 

The main result of our present work is the demonstration of Dissipative Nambu-Hamiltonian 
mechanics of intersecting surfaces as the conceptual framework that underlies strange 
chaotic attractors both in their classical as well as quantum-noncommutative incarnation. 
It reproduces the familiar and well studied attractor dynamics of Lorenz and Rossler in a 
very intuitive manner accounting of their gross topological aspects (double lobe Butterfly 
for the Lorenz system ) or single lobe for the Rossler attractor. 

Quantum Nambu Dynamics of Surfaces 0, El [9], raises also the issue of possible exis- 
tence of Quantum Strange Attractors. Their Quantum behavior was built systematically 
through fuzzifying the classical intersecting surfaces of the ND sector. We demonstrated 
this for the simplest case of the Lorenz system with a linear dissipation. 
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